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ABSTRACT 

The neutrino annihilation is one of the most promising candidates for the jet produc- 
tion process of gamma-ray bursts. Although neutrino interaction rates depend strongly 
on the neutrino spectrum, the estimations of annihilation rate have been done with 
an assumption of the neutrino thermal spectrum based on the presence of the neutri- 
nospheres, in which neutrinos and matter couple strongly. We consider the spectral 
change of neutrinos caused by the scattering by infalling materials and amplification 
of the annihilation rate. We solve the kinetic equation of neutrinos in spherically sym- 
metric background flow and find that neutrinos are successfully accelerated and partly 
form nonthermal spectrum. We find that the accelerated neutrinos can significantly 
enhance the annihilation rate by a factor of ~ 10, depending on the injection optical 
depth. 

Key words: acceleration of particles - accretion, accretion discs - gamma-ray burst: 
general - neutrinos - radiative transfer 



1 INTRODUCTION 

The neutrinos play a very important role in the extremal 
condition in astrophysics. Although photons are strongly 
coupled with the matter in the dense material, neutrinos 
are able to escape due to its weak coupling with the mat- 
ter. Therefore, neutrinos can be an important cooling source 
and significantly affect the dynamics. In addition, they can 
even be a heating source in several cases, e.g., the neutrino 
capture pro cess in the core-collapse supernovae (CCSNe) 
(|Bethelll990t ) and the neutrino annihilation in the central en- 
gine of gamma-ray bursts (GRBs) (|MacFadven fe Wooslevl 
Il999l ). In the hot (temperature T ^ lMeV) and dense (den- 
sity p ^ 10 11 g cm -3 ) gas, even neutrinos are trapped and 
thermalized so that the temperature of neutrinos becomes 
coincident with that of the matter. The surface where the 
neutrinos are decoupled from the matter is called "neutri- 
nosphere" , which is analogous to a photosphere of photons. 

Due to the existence of neutrinosphere, the spectrum 
of neutrinos is often assumed to be a thermal (Fermi-Dirac) 
distribution. However, we know, by the numerous studies of 
photons, that the radiation spectrum can easily be deformed 
from the thermal one in the propagation regime. The non- 
thermal radiation spectrum can be produced by the non- 
thermal spectrum or the different temperature of the scat- 
tering bodies. As for the case of photons, the processes that 
produce the nonthermal spectrum by electron scattering 
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are c alled thermal Compton (see, e.g.. iRvbicki fc Light 



197S ) and bulk Compton processes (jBlandford fc Payne 
19811 ). The former process is driven by the electrons with 
different temperature from photons, while the latter is in- 
duced by electrons with the inhomogeneous velocity field. In 
the case of neutrinos, the bulk motion of scattering materi- 
als could lead to the similar effect and produce nonthermal 
component of neutrinos, which was not investigated so farQ 
Neutrino interactions with matter strongly depend on 
the energy of neutrinos, i.e., the cross section u a with 
£„ being the neutrino energy. Thus, little difference of spec- 
trum (especially at the high-energy region) could lead sig- 
nificantly different dynamics. The neutrino capture and the 
neutrino annihilation are critically important for the shock 
revival of CCSNe and the jet production of GRBs0 respec- 
tively. As for CCSNe, there are significant efforts for solv- 
ing Boltzmann equation of neutrinos with hydrodynamics 



1 Indeed, the term which is related to this effect is included in the 
numerical simulations that solve the neutrino Boltzmann equa- 
tion. However, due to the small radial velocity in the postshock 
region, this effect plays significantly minor role in the context of 
the core-collapse supernovae. Thus, there was no study focusing 
on this effect. On the other hand, we consider the free-fall back- 
ground flow without the shock in this paper so that the radial 
velocity is large enough to make nonthermal component. 

2 Note that there are other alternatives for jet production mech- 
anism of GRBs. Among them, another promising candidate is the 
Poynting dominated jet, which is driven by the magnetic field. 
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numerically because the explosion mechanism (particularly 
delayed-explosion scenario) tremendously relies on neutrino 
physics. Since the neutrinospheres for CCSNe are basically 
spherical (deformation is moderate if the rotation is not very 
rapid), the neutrino transfer can be solved with spherical 
symmetric background, which is reachable even for the cur- 
rent computer resources. In fact, hydrodynamic simulations 
together with neutrino Boltzmann equation have been done 
in spherical symmetri c case (e.g.. iLiebendorfer et aT]|200ll ; 
ISumivoshi et al.ll2005l )FI On the other hand, the central en- 
gine of GRBs is essentially asymmetric (e.g., the compact 
object and the accretion disk system) so that neutrino trans- 
fer should be solved with multi-dimensional treatment. Al- 
though there are a few attempts to solve the neutrino radia- 
tive transfer in multidimensional m anner (e.g.. |Pessart et ail 
l2009l ; ISumivoshi fc Yama da 2012), the long-term dynamical 
simulation is still too computationally expensive so that the 
numerical solutions of full Boltzmann equation are not ac- 
cessible at the moment. Because of these facts, the neutrino 
interactions in the central engine of GRBs are introduced 
with plenty of assumptions. One of them is the thermal spec- 
trum. 

Among the accretion disk models, the neutrino- 
dominated accretion flow (NDAF), in which copious neutri- 
nos are emitted and dominates the cooling, is often discussed 
as a c andidate of the central engine of GRBs. IPopham et alj 
l|l999l ) derived the disk structure and neutrino luminosity 
by solving the group of equations of state, hydrodynam- 
ics, thermodynamics and microphysics in detail. The energy 
conversion efficiencies become extremely high (the annihi- 
lation luminosity L vi , becomes as large as ~ 10 53 erg s _1 ) 
for the mass accretion rate M = WMq s _1 . However, their 
results are too optimistic, as they ignored neutrino opacity 
and ov erestimated the neutrino luminosity. |Pi Matteo et al.l 
|2002l ) showed that the effect of neutrino opacity becomes 
significant for M > IMq s _1 , and recalculated the anni- 
hilation rate including the concept of neutrinosphere. They 
demonstrated that the L vi , increases up to its maximum 
value of ~ 10 50 erg s _1 at M « IMq s _1 and decrease for 
larger M. Thus they concluded that the neutrino annihila- 
tion in NDAF is not a s ufficient mechan i sm fo r liberating 
large amount of energy. iNagataki et al.l |2007l ) performed 
axisymmetric simulation of the collapsar and found that the 
neutrino annihilation is less important than neutrino cap- 
ture as a heating source. These negative results could come 
from the assumption employed in their calculation, i.e., they 
neglected the neutrino emission from the region with neu- 
trino optical depth r Si 1 and employed a thermal distribu- 
tion with a single temperature for the neutrino spectrumQ 



3 Recently, multi-dimensional hydrodynamic simulations with 
neutrino transfer have been per f ormed by several groups 
llBuras et alj |2006| ; iBurrows et"aH |2006| ; iBruenn et all |2009| ; 
ISuwa et al.l hold) by employing several assumptions to re- 
duce the computational costs. More recently, the devclop- 
ment of a full seven-dimens ional Boltzmann solver is reported 
llSumivoshi fc Yamadall2012lV 

4 There are a large number of attempts to amplify the neutrino- 
annihilation rate. Especially, the effects of (general) relativ- 
ity are paid attention, such as beaming b y relativistic motion 
and bending by black- hole space time (e.g. , [Xsano fc Fukuyamal 
120001 , 1200]]) . The relativistic effects for the structure of 



In this paper, multiple scattering of neutrinos and the 
acceleration (i.e., up scattering) in a fluid flow is considered. 
We investigate the impact to the pair-annihilation rate by 
accelerated component of neutrinos. Note that although in 
this paper we consider the parameter regime for the col- 
lapsar scenario that is one of promising candidates of long- 
duration GRBs, the neutrino acceleration process, however, 
is viable for the central engine of short GRBs. Thus, the 
following is applicable both for long and short GRBs. The 
rest of this paper is organized as follows. In iJU W e briefly 
review the concept of neutrinosphere. In >j3l we describe the 
radiative transfer equation for neutrinos and show that its 
solution contains nonthermal component. In SJH we investi- 
gate the effect of nonthermal component of neutrino on the 
neutrino-annihilation rate. We summarize our results and 
discuss their implications in ij5] 



2 NEUTRINOSPHERES 

There are three types of neutrinosphere, w hich are deter - 
mined by the different micro processes (see lRaffeltll200il ). 
Here, we explain these "neutrinospheres" one by one: 

• Number sphere: The optical depth by the emission 
and absorption of neutrinos is about unity. As for v c and v e , 
the electron/positron capture and its inverse process (i.e., 
^ c +n o- p+e~ and £> +p <B> n+c + ) arc important processes. 
As for fx, which represents heavier leptonic neutrinos and 
their antineutrinos (i.e., iv, ^/j> and u T ), the pair pro- 
duction/annihilation processes (i.e., vv -H> 77, vv <-> e + e~, 
NN «-* NNvv) determine the opacity. 

• Energy sphere: The inelastic scattering by electron 
is important process here. The electrons receive the energy 
from neutrinos because the electron rest mass energy (511 
keV) is much smaller than the typical neutrino energy (~ 10 
MeV), which is determined by the matter temperature at the 
number sphere. Inside the energy sphere the neutrinos are 
thermalized due to energy transfer with electrons, which are 
tightly coupled with baryons. 

• Transport sphere: Beyond the energy sphere, the 
elastic scattering by nucleons and nuclei is dominant source 
of the opacity. Because the rest mass energy of these par- 
ticles is much larger than neutrino energy, these scattering 
can be treated as the elastic scatterings 

As for v e and v e , all neutrinospheres provided above are 
almost coincident so that the spectrum is almost thermal. 
On the other hand, vx has distinct radii of neutrinospheres 
|Raffeltll200ll ). Therefore, vx could have nonthermal compo- 
nent, which would be produced between energy and trans- 
port spheres. 



NDAF is also discussed rec ently l lChen fc Beloborodovl 1 20071 : 
IZalamea fc Beloborodovll201ll) . 

5 Note that Raffclt ( 2001) investigated how the recoil term affects 
the spectrum of vx- 
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3 THE RADIATIVE TRANSFER EQUATION 
AND ITS SOLUTION 

Here we consider the neutrino radiative transfer on the back- 
ground of hyperaccreting matter. We focus on the region 
between the energy sphere and the transport sphere, where 
the scattering by nucleons dominates the opacity. 



3.1 Accretion Flow 

Here, we briefly describe profiles of matter as a background 
of neutrino radiative transfer. In this calculation, we employ 
the spherically symmetric accretion flow in order to mimic 
the collapsar. 

Let M be the rate at which matter is accreting, and let 
its radial inward speed be 



u{r) 



1/2 



(1) 



where c is the speed of light and r 3 is the Schwarzschild 
radius. This free-fall velocity profile makes the implicit as- 
sumption that the radiation force on the accreting matter is 
insignificant or equivalently, that the escaping luminosity is 
much less than the Eddington value. The scattering optical 
depth of the flow from a radius r to infinity is given by 



j dr'n(r')a{e v ) = m 



1/2 



C2) 



where n(r) is the nucleon number density (n(r) = 
M /[4irr 2 m p u(r)]), a{e v ) is the scattering cross section for 
neutrino energy e v , and the dimensionless mass accretion 
rate rn = M/MEdd^i respectively. Here, MEdd,„ is the Ed- 
dington accretion rate defined by 



Me 



<j(e u )c 



(3) 



where LEdd,* is the Eddington luminosity of neutrinos, M 



is the mass of the central object, 



is the proton mass, 



and G is the gravitational constant. Since m depends on the 
neutrino energy, we introduce rhkr, which represents m for 
e v — kT with k and T being Boltzmann constant and the 
matter temperature. 

For supercritical accretion into black holes (m 1), it is 
evident from Eq. ([2]) that there should be regions in the flow 
where the neutrinos propagate diffusively. For our problem, 
it is convenient to use not r sc but the effective optical depth 



3 . r s 3 . / £„ \ 2 r s 
_ m _ = _ mfeT (_j _. 



(4) 



This value will replace the radial coordinate r in the ra- 
diative transfer equation in £13.21 It should be noted that 
r depends on e v as well due to the energy dependence of 
the scattering cross section. In addition, we introduce the 
dimensionless neutrino energy as 



kT' 



(5) 



3.2 Radiative Transfer Equation 

Now we proceed to write down and solve the radiative trans- 
fer equation. For the kinetic equation, we start from equa- 



tion (18) of iBlandford fc Payne! { 198lh for the neutrino oc- 
cupation number f v (r, e„), 

3/t(e„) J ( 6 ) 



dt 



+ u ■ V/, 



+±<v 



'dT u + ^ e * 



where «(r) = n(v)a{e v ) is the inverse of the scattering mean 
free path and j(r, e„) is the emissivity. Here, we neglect the 
recoil term, which affects the neutrino spectrum only for 
the regime of e v £ m p c 2 ~ 1 GeV (typical neutrino energy 
is ~ 10 MeV). Substituting the inflow velocity u = — «(r)f, 
where f is the radial unit vector, and taking into account 
the spherical symmetry, Eq. ([6]) becomes 



dfv 

dt 



d 2 U 



- 2r 



df» 
8t 



1 dfy 

2 X dx 



+ 



3ck(x) ' 



(7) 



where t = 3cnt/r is the dimensio nless time. The dimension - 
less spectral energy flux F{r,e v ) (|Blandford fc PavneHl98lh 
is written in the new variables as 



F(t, x) oc x 



2t 
3m H 



1/2 



2t 



+ 1 



df 



-+ l -x d -± 

dr 3 dx 



(8) 



Note that we are interested in the spectrum at a certain 
radius, not the optical depth, which depends on both the 
radius and neutrino energy. Thus, by combining Eqs. Q 
and ©, we evaluate the spectral energy flux at a certain 
radius as 



F(t,x) 



2_L + 1 )dfv 
3 J dr 



+ 1 -x^ 

3 dx 



(9) 



3.3 Analytic Solution 

In this subsection, we neglect the last term in the right hand 
side of Eq. in order to obtain an analytic solution. This 
is because this term changes the number of neutrinos, which 
would have minor contribution betwee n the energy sphere 
and tr ansport sphere (see ij2}. Following lPavne fc Blandfordl 
(|l98ll ), we solve Eq. Q using variable separation with the 
form 



Ut,x) = R{t)t 



5/2, 



(10) 



Here a becomes an eigenvalue of the following confluent hy- 
pergeometric differential equation 



d 2 R 



2 J dr 



a - 10 



R = 0. 



(11) 



The physical solution of Eq. (|1 1 f) fulfills a constant spectral 
flux of neutrinos as t -> and adiabatic compression of 
the neutrinos r — > oo. The relevant solution can be evalu- 
ated as an infinite sum of generalized Laguerre polynomials 
L„ (2t). The corresponding eigenvalues a n are given by 



a„ = 4n + 10; n = 0, 1, 2 
These values are different from I Payne 



(12) 

Blandfordl (|l98lh 



because the cross section depends on energy for the cur- 
rent case. Note that this spectral index does not imply the 
observable spectrum at a certain radius because r depends 
on not only a radius but also the considered neutrino energy 
(see Eq.|4|). In order to obtain a spectral flux at r, we should 
multiply Eq. (£9j) by x 3 because df v /dr\ r oc r 3/2 oc x 3 r~ 3/2 
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for t — > (the other terms drop faster than this term). Thus, 
the hardest spectral component (i.e., n = 0) of F becomes 
oc e~ 4 . 

In principle, the global solution of Eq. (0 can be given 
by summing up infinite series expressed by L% 2 (2r)x~ arl 
with coefficients determined by the boundary condition. In 
fact, iPavne fc Blandfordl (|l98ll ) gave the analytic solution 
with the delta-function distribution function at the injec- 
tion optical depth, in which all coefficients are expressed 
by generalized Laguerre Polynomials and Gamma functions 
(see Eqs. 10 and 11 in their paper). However, the analytic ex- 
pression for the arbitrary boundary condition is not always 
representable using known functions. Thus, in the following 
we solve Eq. (Q) numerically with the thermal distributions 
of neutrinos at the energy sphere as a boundary condition. 



Spectrum w/ bulk term — 
Spectrum w/o bulk term 




x 



3.4 Numerical Solution 

In this subsection, we present our numerical solution of Eq. 
(Q) . The last term is omitted again because the interested re- 
gion is between the energy sphere and transport sphere (see 
0. We use the relaxat ion method for the boundary problem 
(e.g., iPress et al.lll992T ). in which the stationary solution is 
achieved by infinitely long exposure of the time dependent 
equation. 

In solving Eq. 0, we change this equation to a finite- 
difference form using 



d 2 U _ 2 
dr 2 Srt + ST i+1 

dr 



fi,3 fi,j _ f 
J v J v J I 



i— 1,J 



5n 



113) 
(14) 

2~ J Sxj+dxi+i ' (15) 

where i and j denote the grid point of r and x, respectively. 
The grid points are determined by the rule as 



1 df» 
2 X dx 



Sri + 5r i+ i 



f»..j+i _ ft.: 

.IV .IV 



n = Ti-i + 5n, 

Xj = Xj-i + Sxj, 
5n = r T Sn-i, 
Sxj = r^Sxi-i, 



(16) 
(17) 
(18) 
(19) 



where ry and r x are constants larger than unity. We set 
Sti/tl — Sxi/xi = 0.02. The calculations are performed on 
a grid of 200 zones for r from 0.01 up to t and 500 zones for 
x from 0.1 to 100. A test calculation and comparison with 
an exact solution are given in Appendix. 

The boundary condition is given at t as 



fv{r ,x) = 



^5/2 



(20) 



x 5 1 + e x 

where the factor Tq^ 2 /x 5 means the correction, which leads 
to the thermal distribution function at a radius r. In this 
study, we set the above boundary condition with only one 
parameter t for simplicity. In order to make more realistic 
boundary condition, we should consider the microphysical 
processes that change the neutrino number and energy in 
detail, which is beyond the scope of this paper. Here we do 
not care about the normalization factor because all equa- 
tions solved in this study are linear to /„ as we omit the 
source term j. Needless to say, we should care about the 



Figure 1. The emergent spectral energy flux estimated by the nu- 
merical solution. The boundary condition is given at to = 5. The 
correction for the conversion from r-space to r-space is included 
(see text for details). The red solid line is the solution of the full 
equation of Eq. J7), while the black dashed line is the solution 
of the kinetic equation without the bulk term (i.e., thermal spec- 
trum). The grey dotted line represents the power-law spectrum 
of £^ 4 , which is the analytic solution (see text for details). 



normalization with detailed source term because neutrino is 
fermion so that there is a significant effect by Pauli blocking 
for / ~ 1. 

In Figure [1] we show the numerical solution of the di- 
mensionless spectral energy flux obtained by solving Eqs. 
and (j9} ■ The red solid line represents the emergent spectrum 
of full equation and the black thin dashed line is spectrum 
obtained by the kinetic equation without bulk term (i.e., 
^(V-u)£v^- in Eq.|6]), that is, a thermal spectrum. One can 
see that neutrinos are upscattered by the infalling material 
and the nonthermal spectrum is generated. As indicated by 
grey-dotted line, the emergent spectrum is power law with 
e~ 4 for x = Ev/kT ^ 10, which is consistent with the ana- 
lytic solution with n = obtained in the previous section. In 
this calculation, the boundary condition is given at tq — 5, 
where /„ has a thermal distribution with a temperature, T|j 
The flux is estimated at r = 0.01, where the spectral evolu- 
tion is almost completed. The normalization of both spectra 
is determined by the total number flux, J(F/x)dx, being 
unity. 

We show the different solutions with different To in Fig- 
ure [2] It is obvious that higher to leads to harder spec- 
trum. It should be noted that the position of the energy 
sphere depends on the elementary process such as ve «-> ue, 
e + e~ -f-> vv, and NN bremsstrahlung. The first process is 
related to the thermalization and the others are related to 



e According to numerical simulations of core-collapse super- 
novae, which include detailed microphysics and the radiative 
tran sfer, the temperature of vx ranges from 4 MeV to 10 MeV 
(see iHoriuchi et al.l |2009| . for a collective reference of recent nu- 
merical simulations). Thes e values can be used in the case of long 
GRBs. As for short GRBs. lSetiawan et all j2006T) showed that the 
average energy of vx range from ~ 5 MeV to ~ 27 MeV, corre- 
sponding to the temperature from ~ 2 MeV to ~ 9 MeV with 
vanishing chemical potential, similar to values of CCSNe. 
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Table 1. Properties of numerical solutions 
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Figure 2. The emergent spectral energy flux with different to 
indicated by different lines. The higher To leads to the harder 
spectrum due to efficient up-scattering by the infalling material. 



both the thermalization and emission/absorption. The dom- 
inant thermalization process depends on the background 
fluid temperature, density, and abundance, which are much 
beyond the scope of this paper. Thus, we simply parametrize 
the injection tq and see the dependences on it (see also 



the injection 
lRaffeIul200lh 



4 NON-THERMAL NEUTRINOS AND THEIR 
ANNIHILATION 

Now we move on to estimate the neutrino annihilation rate, 
which strongly depends on neutrino energy. The energy de- 
position r ate via neutrino-annihilation (v + v — > e + + e~) is 
given by (|Goodman et al.lll987l ; ISetiawan et al.ll2006t ) 



E v o — CF3 „F3 n 



(el)(eu) + (ef)(e v ) 



(21) 



where F i>v = / /„ej,de„, (e„) = F 3?I/ /F^ v , and (el) = 
Fi tV /F2,ui respectively. The factor C includes the weak in- 
teraction coefficients and the information of the angular dis- 
tribution of the neutrinos so that to calculate this factor 
we should determine the geometry of the neutrino-emitting 
source. Since this factor is expected not to change signif- 
icantly by including the neutrino acceleration process, we 
concentrate on the effect of the spectral change from here. 
For simplicity, we assume the spectrum of v and v are iden- 
tical. Then, we get 



E v „ oc 



FL(el) 
(e v ) ■ 



(22) 



We can evaluate the amplification of the neutrino annihila- 
tion rate by the accelerated component of neutrino produced 
by the bulk motion of background matter using F3,„, (e v ), 
and (el). By assuming that the neutrino number flux (-Fb.i/) 
does not change by including this effect, we get Fs tV oc (e v ), 
then E V v oc (e„)(el) Therefore, we can evaluate the ampli- 
fication only by (ej) and (el). 

Table [T] shows the integrated values of emergent spec- 
trum. It is obvious that both the mean energy (e v ) and the 



(e^ 



(4) 



At 



■A-the 



0.1 


200 


500 


1.01 


1.02 


1.03 


0.2 


200 


500 


1.03 


1.05 


1.08 


0.5 


200 


500 


1.07 


1.16 


1.24 


1.0 


200 


500 


1.16 


1.37 


1.59 


1.5 


200 


500 


1.26 


1.65 


2.08 


2.0 


200 


500 


1.37 


1.99 


2.73 


3.0 


200 


500 


1.60 


2.83 


4.52 


5.0 


200 


500 


1.95 


4.49 


8.77 


7.0 


200 


500 


2.18 


5.72 


12.5 


10.0 


200 


500 


2.43 


7.12 


17.3 


1.0 


500 


1000 


1.16 


1.36 


1.57 


1.0 


750 


1500 


1.16 


1.37 


1.59 



t Numerical grids for t. 
' Numerical grids for x. 

§ The average energy of emerged spectrum. 

^ The amplification of the neutrino-annihilation rate (Eq. 122 I t 

compared to the thermal distribution. 



mean-square energy (el) increase compared to the thermals 
spectrum due to up-scattering by the infalling materials. As 
a result, the neutrino annihilation rate is significantly am- 
plified by the accelerated component. In addition, we show 
the convergence check with higher resolution in this table 
(see last two lines). Due to much more expensive numerical 
cost, we just calculate the model with to = 1 and confirm 
the validity of the lower resolution calculation. 



5 SUMMARY AND DISCUSSION 

In this paper, we consider the spectral change of neutri- 
nos induced by the scattering of the infalling materials and 
amplification of the annihilation rate, which is one of the 
well-discussed jet production mechanism of GRBs. We solve 
the kinetic equation of neutrinos in spherically symmetric 
background flow and find that neutrinos are successfully ac- 
celerated and partly form nonthermal spectrum. We find 
that the accelerated neutrinos can significantly enhance the 
annihilation rate by a factor of ~ 10, depending on the in- 
jection optical depth. 

In this study, we tried to demonstrate the effect of the 
up-scattering by the bulk motion of the material and just 
assumed the injection optical depth with parametric man- 
ner. More realistic injection is obtained by the insight of the 
energy sphere, whose position is determined by the neutrino- 
electron inelastic scattering as follows. The scattering opac- 
ity at energy sphere, where the inelastic scattering with elec- 
trons freeze out, can be calculated by the ratio of cross sec- 
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tio: 



: (r es ) 



2e„ 



n e iTe 3Y c fcT(r es ) 



(23) 



where r ca is the optical depth of electron inelastic scattering, 
r cs is the radius of energy sphere, «n and n c are the number 
density of nucleons (neutrons and protons) and electrons, 
and Y e = n e /n^ is the electron fraction, respectively. Since 
r es is definitely 2/3 at the energy sphere, 

4e„ 



Z (^*es) 



9Y B kT(r c 



(24) 



The typical temperature of neutrinospheres is ~ 4 MeV 
ljjankall200ll ) so that r ac is larger than 2/3 for neutrinos of 
the energy e v ^ 0.6(Y C /0.1) MeV. The typical value of Y e is 
~ 0.1 at the region where the electron capture is significant 
so that almost all the neutrinos are trapped by the nucleon 
elastic scattering at the energy sphere. Although the injec- 
tion optical depth in this study should be energy dependent 
as shown above, we neglect this effect for simplicity. 

Next, we discuss about the neutrino species. In this pa- 
per, we consider the region between the energy sphere and 
the transport sphere, i.e. the optical depth is larger than 
unity for the neutrino-nucleon elastic scattering. Note that 
usually these neutrino spheres are coincident for v c and v e 
due to the presence of the charged current for these neutri- 
nos so that the neutrino acceleration studied in this paper 
is possible only for v^, v T and their anti particles. How- 
ever, for the case of neutron number density being much 
larger than protons' one, the charged current reaction of v e 
(p e +p — > n + e + ) is negligible so that the reactions relevant 
to i> e become similar to those of heavier leptonic neutrinos. 
The transport opaci ty for the ne utral current scattering pro- 
cesses are given by Jjankal l200lll 



5a 2 + 1 ctq(4) p 
24 (m c c 2 ) 2 m u 



(Yn + Y p ) 



(25) 



Here m u ~ 1.66 x 10~ 24 g is the atomic mass unit, a = 
— 1.26, and Y n = n n /njv and Yp = n p /riN are the num- 
ber fractions of free neutrons and protons, i.e., their particle 
densities normalized to the number density of nucleons, re- 
spectively. In cases of v c and v c also the charged-current 
absorption on neutrons and protons, respectively, need to 



7 The total cross section of neu trino-electron inelastic scattering 
is lBurrows fc Thompson! i2002l') 

3 e v kT 
-o-Q- 



" o u I 212 ' 

where <ro ~ 1.7 X 10 — 44 cm 2 is reference neutrino cross section, 
m c is the electron mass, T is the temperature of electrons. On 
the other hand, the total cross section of neutron-neutrino clastic 
scattering is 

/ \ 2 

°0 



The reason why we employ the cross section for neutrons is that 
due to the electron capture (p + e - — > n + u e ) the neutron frac- 
tion increases, whereas the proton fraction decreases inside the 
ncutrinosphere. Therefore, neutrons are dominant target parti- 
cle for propagating neutrinos. However, it should be noted that 
the total cross section of proton-neutrino scattering differs from 
neutron only ~ 20%. 



be taken into account due to their large cross sections. The 
absorption opacity is (Janka 2001) 



3a 2 + 1 cr (e 2 ) 
4 (m c c 2 ) 2 



Yn 

Yr> 



(26) 



From Eqs. (|25|) and (|26[) the scattering dominates the opac- 
ity for v e provided Y p < .26. In this case, the tra nsport 
sphere and number sphere (Janka 1995; Raffelt 2001) sepa- 
rate from each other for v e . The time stationary solutions of 
hyperaccreting flow imply tha t Y e can be as small as ~ 0.1 
|Kawanaka fc Mineshigei 20071 ) . similar to the case of core- 
collapse supernova (just above the neutrino sphere, Y e ~ 
0.05 - 0.1). Therefore, it is expected that the acceleration of 
v a would naturally occur in the collapsar system. 

Our finding suggests that the detectability of MeV neu- 
trinos is also enhanced because the expected detection num- 
ber oc -f2,^(e 2 )- If this neutrino acceleration works only for 
ux , the neutrino oscillation would produce v a , which is main 
observ able for waiter Cerenkov detectors. ISuwa fc Murasel 
(|2009l ) estimated the expected number from the hyperac- 
creting accretion flow with thermal spectrum of kT = 3 MeV 
and argued that GRBs are observable for £ a few Mpc by 
Super-Kamiokande and several Mpc by Mton detector. The 
neutrino acceleration process thought in this paper would 
push out the detectable horizon of MeV neutrino farther. 

At last, we comment on our assumptions in this study. 
Firstly, we employed the diffusion limit for whole region, 
which is not valid for optically thin rsgion, Tsc < 1. The 
conclusion, however, does not change if we somehow include 
effects of optically thinness, since the spectral evolution is 



determined by the region t s , 



1. Secondly, we dropped 



out the recoil term from the kinetic equation. The average 
energy of neutrinos is not affected by this because the re- 
coil term changes the spectrum for ^ 1 GeV which is much 
higher than the typical energy of neutrino spectrum. Thus, 
if we include the recoil term (that is much more complicated 
than terms included in this study), the conclusion does not 
change very much. Thirdly, we fixed the background matter 
flow as free fall. The energy gain of neutrinos must come 
from the matter so that the back reaction should be in- 
cluded when the total energy of neutrinos reaches as large 
as the matter kinetic energy. To compare these two quan- 
tities, more detailed source term is necessary. The final an- 
swer can be obtained by solving neutrino radiation hydro- 
dynamic equations in self-consistent way, which is far be- 
yond the scope of this simple study. Fourthly, we omitted 
the relativistic effects, e.g., the Doppler shift and gravita- 
tional redshift. In order to include these effects, we should 
reformulate by covariant formulation of the radiative trans- 
fer i^sdf^con^tentway. One can find such a formulation 
in IShibata et all (|201ll ). Finally, the spherical symmetry is 
also one of the largest assumption in this study. This as- 
sumption is partially valid because NDAF solution has large 
disk height due to the large contribution of gas for pressure 
(otherwise if the disk is supported by the centrifugal force, 
the disk height is almost negligible compared to the disk ra- 
dius). Considering the disk structure, neutrinos can escape 
from the accretion flow to the vertical direction before the 
spectral evolution completes. Whether this effect amplifies 
or suppresses the neutrino annihilation is not trivial because 
there are two possible opposite effects; neutrinos emitted at 
deep position, which are experienced significant acceleration, 
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Figure Al. Absolute value of the relative error of the numerical 
solution in the case of tq = 1. 



can more easily escape than the spherically symmetric con- 
figuration, while the acceleration might not complete due 
to the earlier escape. To give more concrete result, more 
detailed calculation is strongly required. 
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APPENDIX A: TEST CALCULATION 

In this section, we show a test result of our numerical cal- 
culation. A special solution of Eq. © is given by 



r , \ 5/2 - 

f v (r, x) — T x 



(Al) 



with j = and adequate boundary conditions. In Fig. 
IA1I we show the absolute value of the relative error, |1- 
(numerical solution) /(exact solution)] in the case where 
ro = 1 as an example. The grid setup is the same as one 
used in Sec. 13.41 We find that the relative error is always 
less than ~ 10~ 2 for (IA1|) . Hence, the error does not affect 
the conclusion in ij4] 
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